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A  TRANSPORTATION  PROBLEM  INVOLVING 
SOURCE-LOCATION  OPTIMIZATION 

R,  L,  Slelken  Jr, 

A  generalisation  of  the  capacitated  plant-location  problem  is  formulated 
and  several  possible  modifications  noted.  Four  solution  techniques  are 
briefly  discussed:  exhaustive  enumeration,  probabilistic  search,  zero-one 
mixed  integer  linear  programing,  and  an  Iterative  procedure.  The  relative 
attractiveness  of  the  iterative  procedure  is  illustrated  in  three  examples. * 


1,  Introduction 


The  general  source-location  problem  under  consideration  is  how  to  supply 
J  destinations  with  D^»  B£»  . ..,  units  from  K  possible  sources  at  a  minimum 
cost  when  the  K  sources  have  capacities  B^,  and  any  subset  of  the 

sources  can  be  located  at  any  one  of  I  locations.  A  minimal  cost  solution 
involves  the  specification  of  each  source's  location  and  the  allocation  of  the 
demands  D^»  ,  Dj  among  the  sources. 

A  aero-one  mixed  integer  linear  programing  formulation  of  this  problem 
is: 

^  1  ^*1  Ei*l  Ej«l  ckij  *kij  +  kE1  J.  \i  V 
gsFject  to  the  constraints 

Ek»i  ELl  *kij  *  Dj*  j  "  1 . J’ 

Ej=l  *kij  -  “ki  V  1  "  *•  •**  1  md  k  "  U  K* 

Ei®l  \i  —  1*  k  “  1*  *•**  K» 

u^  «  0  or  1,  i  «  1,  ...»  I  and  k  ■  1,  ...»  K, 

^  Oj  all  k»  i,  j» 

where 

.1  if  source  k  is  at  location  i, 
uki  "  ''O  otherwise; 
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fki  a  the  fixed  cost  for  source  k  being  at  location  i; 
x^  -  the  number  of  units  transported  by  source  k  from  location  i  to 
destination  j;  and 

c,^  «  the  cost  of  source  k  producing  one  unit  at  location  i  and  transporting 
it  to  destination  j. 

The  equality  in  (2)  implies  that  the  demand  at  each  destination  is  satisfied. 

The  constraint  in  (3)  insures  that,  if  source  k  is  at  location  i,  then  the  capacity 
of  source  k  is  not  exceeded,  and,  if  source  k  is  not  at  location  i,  then  nothing 
is  transported  by  source  k  from  location  i.  The  inequality  in  (4)  implies  that 
each  source  it.  located  in  at  most  one  location.  The  constraints  also  imply  that, 
if  a  source  is  no£  located  anywhere  all,  then  it  cannot  help  satisfy  the  demands. 

Of  course,  the  physical  interpretation  of  the  sources,  destinations,  costs, 
etc,  1b  quite  unrestricted.  For  example,  the  sources  could  be  refineries,  the 
dsotinations  could  be  military  installations,  and  the  coats  in  terms  of  time. 
Alternatively,  the  sources  could  ba  generators  of  electricity,  the  destinations 
cities,  and  the  cobcs  monetary. 

2.  Scats  Possible  Modifications  of  the  General  Source-Location  Problem. 

Restrictions  concerning  the  feasibility  of  certain  source-location  combinations 
can  be  easily  incorporated.  For  example,  if  at  most  sources  can  be  located  at 
location  i,  then  the  problem  should  be  modified  by  adding  the  constraint 

K 

z  v  <  l  .  m 

k-1  1  4 

Of  course,  if  source  k  can  only  be  located  at  a  location  in  a  subset  of  the 
I  locations,  than  the  problem  should  be  modified  by  deleting  all  variables 


V 


-3~ 


and  x^  with  i  |  T^ 

The  extension  of  the  probleta  to  a  situation  involving  more  thtn  one  product 
is  conceptually  si  mole.  If  there  are  P  products  ad  a  subscript  p  is  used  to 
inc. -cate  the  p-th  product,  the  extended  problem  is 


4 

■nin  . 

y  i  J  N  i 

l  ,  £j  .  Lj  .  c~4 .  X.  +  E  Eu,.  f.,, 

jc»?  i-1  j“l  rfijp  TujP  Ual  it>1  Ki  ki 

(8) 

subjec  to  thf 

■r‘  Mvraints 

r* 

k“l 

Ei»l  \ijp  "  Djp*  allj»P’ 

(9) 

l3 

>1 

*ki;)P  -  “ki  Bkp»  a11  k»  p» 

(10) 

E1 

£i«l 

“ki  Ji  1,  k  ■  1 . K, 

(11) 

v  ■  0  01  l>  811  k*  lf 

(12) 

allk>  “• 


(13) 


Of  course*  the  general  source-loeation  probiea  as  formulated  in  (l)-(6) 
encompasses  the  special  case  ia.  which  all  K  source*  are  identical*  However,  if 
all  K  sources  are  identical,  the  problem  coo  else  be  formulated  ae 


Ein  Zi-1  £j-l  Cij  2ii  +  J*  Vi£i 


(U) 


•4> 


.  object  to  the  constraints 


1, 


• « •  ^ 


€15) 


B, 


i  * 


^  ; 


i» 


Ci6> 


V£  IK, 


cm 


“  0|  If  2,  ■  •  <  t  Kj  1  "  1*  o » •  >  X*  (Is) 

Xy  >  0,  all  i,  j,  (19) 


share 

B  *  the  capacity  of  a  source: 


*  the  cost  of  producing  one  unit  at  location 


destination  j; 


and  transporting  it  to 


*  the  number  of  units  produced  at  location  i  and  transported  to 
destination  j; 

-■»-  the  fixed  coat  per  source  located  at  location  i;  and 


*  the  nia&fef  of  sources  at  location  i. 

This  alternative  fotsulstion  suggests  that  for  this  particular  case  specialised 
solution  techniques  should  be  used.  The  deyelopssnt  of  such  specialised  techniques 
ia  in  progress  at  this  ti*s«  - 


-5- 


3,  A  Survey  or  the  Literature 

The  "capacitated  plant-location  problem"  considered  by  Bulfin  and  Unger 
[3],  Davis  [9J,  Ellwein  [12],  Gray  [14],  [IS],  and  [16],  Marks  [19],  Sa  [21], 
apd  Spielberg  [24]  can  be  formulated  as 

T  K& 

»in  \al  zjol  ckj  +  zk»!  fk\ 

subject  to  the  constraints 


Ek-1  *kj  "  Dj*  j  “  X»  J 


*kj  -  \  V  k  *  1 . K* 


-  0  or  1,  k  *  1,  K* 


L> : 


>  0,  all  k,  j 


Thus,  from  one  viewpoint,  the  capacitated  plant-location  problem  is  a  general 
source-location  problem  with  I  -  1.  Alternatively,  from  the  opposite  viewpoint, 
the  general  source-location  problem  is  a  capacitated  plant-location  problem  with 
K  sets  of  I  mutually  exclusive  plants.  Some  discussion  of  mutually  exclusive 
plants  is  given  by  Marks  [19]. 

The  "simple  plant-locaticn  problem"  considered  by  Celebiler  [4],  Efroymson 
and  Ray  [11],  Feldman,  Lehrer,  and  Ray  U3]f  Manne  [IS],  and  Spielberg  [22]  and 
[23]  has  the  same  formulation  as  the  capacitated  plant-location,  problem  except 
that  the  capacity  of  each  source  is  assumed  to  be  at  least  equal  to  the  total 
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demsftd,  This  assumption  si'iplifies  the  problem  considerab3y  since  it  implies  that 

each  demand  D4  can  always  be  optimally  supplied  from  the  "nearest"  source. 

J 

Warehouse  location  problems  which  can  be  simplified  to  general  source- 
location  problems  were  considered  by  Balinski  and  Mills  [1] ,  Baumol  and  Wolfe 
[2],  Kuehn  and  Hamburger  [173,  and  Marks  [19]. 

Chapelle  [a],  Cooper  [6],  [7],  and  [8],  and  Raincaek  and  Hartley  [20] 
investigated  the  following  s»r:rce-location  problem:  If 

D_j  =  the  demand  for  a  product  at  destination  j,  j  -  1,  ...»  J; 

(a^b^)  =  the  coordinates  in  two-dimensional  Euclidean  space  of  destination  j; 

=  the  capacity  of  source  k,  k  =  1,  ...,  K; 

(a.  jb,^)  =  the  coordinates  in  wo^dimensional  Euclidean  space  of  source  k; 


* »  the  number  of  units  transported  from  source  k  to  destination  j  ; 
than  the  objective  is  to 


Sk»l  Ej«l  *kj  [UDj  '  V  +  (bDj  **  bk)  3 


2,1/2 


subject  to  the  constraints 


W  \j  *  Dj 


^D. ,  ja‘l*...,J» 


EJ-1  *kj  -  V  k  '  1 . *’ 


bki°-  “kl  -°-  *u  k-  5' 
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This  problem  is  the  same  as  the  general  source-location  problem  except  that 
(i)  the  location  of  each  of  the  K  rces  is  not  restricted  to  one  of  a 
finite  number  of  possibilities,  and 
(ii)  c,  . ,  the  cost  of  source  k  producing  one  unit  at  location  i  and 

transporting  it  to  destination  j,  has  been  replaced  by  the  Euclidean 
distance  from  the  location  of  sourer  k  to  destination  j. 

4,  Some  Solution  Techniques 
4»1  Exhaustive  Enumeration 

To  solve  the  general  source-location  problem  by  exhaustive  enumeration,  the 

optimal  allocation  of  the  J  demandsamong  the  K  sources  would  have  to  be  determined 
£ 

for  each  of  the  I  possible  assignments  of  the  K  sources  to  the  1  locations.  For 
t. -h  assignment  the  determination  of  the  demand  allocation  is  a  capacitated 
plant-location  problem;  that  is 

”ln  4.1  4-1  ckj  *kj  +  4.1  £k  “k  (29) 

subject  to  the  constraints 


Vl  *kj  “  Dj 


D, ,  j51!, 


*kj  -  Bk  V  k  =  1#  ”• 5  K’ 


Ufc  “  0  or  1*  k  "  1 . K 


xRj  >  °,  all  k,  j. 


B?anci>and-bound  treatments  of  this  problem  have  been  given  in  [3],  [9], 
[12],  [21]  sad  [24],  The  empirical  evidence  suggests  that  the  branch-and-b^md 
method  would  be- most  useful  It  small  problems,  around  25  integer  variables  or 
less. 


Thera  are  at  least  two  difficulties  with  the  exhaustive  enumeration  approach; 

K 

(i)  the  number  of  possible  assignments,  I  ,  increases  dramatically  as  I  and 
K  increase,  and 

(ii)  the  only  available  "fast"  algorithms  for  solving  medium  or  large 

capacitated  plant-location  problems  are  approximate  routines  -  see, 
for  example,  [21]. 

4.2  Probabilistic  Search 

Since  the  number  of  possible  source-location  assignments  may  be  quite  large 
in  a  realistic  situation,  it  may  only  be  feasible  to  explicitly  evaluate  a  subset 
of  them.  A  simple  probabilistic  search  procedure  which  explicitly  considers 
only  a  subset  of  the  possible  source-location  assignments  is  as  follows.  Select 
a  random  sample  of  size  n  without  replacement  from  the  set  of  all  possible 
assignments.  Solve  the  n  problem*?  (29}-(53)  corresponding  to  n  assignments 
selected.  Let  z*,  z£,  ...,  s*  represent  the  optimal  values  of  the  objective 


functions  in  these  problems.  Then,  if  z^,  s ^ . z  R  are  the  order  values 

of  ^2 ,  • i  s  ^  with 


*<i)  -  z(i+l) *  i 


1. 


IK-1, 


(34) 


it  follows  that 


1  <  m  <  n 


<  zf  .)  >  1 
a  -  (r)  - 


TK_r  1 

“  >/C 


-  C^rrX&g-l)  ***  (IK-r-n-M) 


IK(IK-1)  '*•  (IK-n+I) 


Thus,  for  any  r,  the  probability  of  obtaining  an  assignment  at  least  as  good  as 


the  r-th  best  assignment  can  be  made  as  close  to  one  as  desired  by  taking  the  sample 


size  n  sufficiently  large. 


One  procedure  for  sequentially  deterging  which  source-location  assigoements 


are  to  be  evaluated  is  as  follows.  Let  a  be  a  constant  such  that  0  <  a  <  1.  Let 


z^,  ...»  *  R  represent  the  optimal  values  of  the  objective  factions  in  the  I 

^  K 

capacitated  plane-location  problems  (29)-(33)  corresponding  to  the  I  possible 


assignments.  Assume  that  a  histogram  of  z^,  ...,  z  ^  can  be  closely  approximated 


by  a  density  function  of  a  given  form.  Take  a  random  sample  of  size  one  from 


the  set  of  I  assignments.  Solve  the  problem  in  (29)-(33>  corresponding  to  the 


assignment  selected.  Estimate  the  parameters  of  the  approximating  density 


function.  Then  repeat  -hese  three  steps  until  a  selected  assignment  makes 


(29)  less  than  the  a-th  percentile  of  the  distribution  function  corresponding  to 


the  estimated  density.  This  procedure  was  discussed  by  Cooper  [7j  in  relation 


to  a  similar  source-location  problem. 


In  comparing  these  two  probabilistic  search  procedures,  the  former  procedure 


has  the  advantage  of  selecting  a  fixed  number  of  source  location  assignments  for 


consideration  and  does  not  necessitate  any  assumption  about  the  distribution  of 


Sp  z^*  *  g*  On  the  other  hmd,  Cooper's  sequential  procedure  could  result 


in  a  smalles?  sample  size. 


■ tm 


1 
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4.3  Zero-One  Mixed  Integer  Linear  Programing 

The  general  source-location  problem  is  a  zero-one  sixes  integer  lines? 
programming  problem:  hence,  any  general  mixed  integer  linear  progressing  algorithm 
may  theoretically  be  used  to  solve  it.  In  particular,  the  relatively  new  hrench- 
and-bound  methods  and  implicit  enumeration  methods  seem  applicable e  However, 
computational  experience  on  the  simple?  capacitated  plant "location  problems 
indicates  that  these  methods  may  only  be  coi^jutatiiraally  feasible  for  relatively 
small  Tobieias. 

4.4  An  Iterative  Algorithm 

If  the  source-location  variables,  u^,  were  known ,  then  the  general 
source-location  problem  would  reduce  to 


,  rK  .T  „I 
sin  E.  -  >.  -  e- 


'kaI  ^j«l  “i«I  ^ij  ^kij 


(36) 


subject  to  the  constraints 


^k*!  £i«l  \ij 


&j  *  j  “  1»  • . • »  J  +  1^ 


(37) 


where 


Z1»l  *kij  "  “ki  Bk’  aU  *•  *■* 

t  °»  a3Ll  k*  U 


B 


E,K  ,  I1 


’j+l  ’  ^k-l  ^i"!  “ki  Bk  “  Ej»l  Dj 


(33) 

(39) 

(40) 


This  problem  is  just  a  simple  transportation  problem. 

On  shs  other  aaad,  if  the  allocations  of  tbs  demands  D^,  ,  .,*  £*j  aasoag  the 
X  sources  wars  known  and 

8^  *  the  number  of  unite  to  be  transported  from  source  k  to  destination  j, 
thso  fth®  general  source- location  problem  would  reduce  to 


'dzm 


-• 


wmss^M 


mm 

-Cv'\'Xcf-§| 

pjlrf 

WM 


x  i 

min  £  £ 

k*l  i«I 


\i  ffki  +  ^  Ckij  °kj* 


sitsbjcct:  to  the  constraints 


e|h1  u^  »  1,  for  all  k  such  that  E^^  >  0, 


£.N 


5T-i!-sJtea®S: 


s-v*; 


ft! 


mmmm 


yj.  •-••-- v  .  o*<J 


4 


»  0  cr  1,  all  k,  i. 


An  optimal  sol-tion  to  this  problem  is  simply 


uJti  »  1  if  i  »  i*(k)  and  E  x  >  0 


“  0  otherwise 

where,  for  each  k,  i*(k)  is  the  smallest  positive  integer  i*  such  that 


^h-4  *  *  ^ 


°kj  "  i  <mf<  j  Uki  +  *  Ckij  °kj]  * 


The  simplicity  of  both  optimally  allocating  demands  among  sources  for  a  fixed 
source- location  configuration  and  determining  an  optimal  source-location  configuration 
for  a  fixed  allocation  of  demands  among  sources  suggests  the  following  iterative 
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schema  for  obtaining  approximate  solutions  to  she  general  source-location 
problem: 

(i)  Chooea  an  initial  source- location  configuration. 

(ii)  For  this  source-location  configuration*  determine  the  optiir  l  allocation 
of  demands  among  the  sources. 

(iii)  For  this  demand  allocation,  determine  the  optisa^i  aourae-iocaticn 

configuration.  Return  to  (ii)  with  this  a**?  configuration. 

This  scheme  of  alternately  evaluating  optimal  demand  allocations  sad  optimal 

source- location  configurations  can  be  formalised  into  the  following  procedure: 

1.  Select  an  initial  value  for  each  u,  .  with 

si 


■  0  or  1,  all  k,  i 


and 


EAi  -  l>  k  •  l> 
i»l 


V 

»  •  •  f  • 


2. 


Sclve  the  transportation  problem  in  (36)-(40)  with  the  u^’s  biased 

their  current  values.  Represent  the  optimal  demand  allocation  thus 

obtained  by  { D.  ;  k  -  1,  ...,  K  and  j  **  1,  ....  J}  whare  D.  ,  is  fb^  nigsbxr 
kj  KJ 

of  units  to  be  transported  from  source  k  to  destination  Jj  s.e? 


I 

E 


'  A  %u- 


Generate  new  values  in  accordance  with  (44),  so  that  the  corresponding 
source-location  configuration  is  optimal  for  the  new  allocation 

determined  in  stop  2. 


3. 
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4.  If  the  aoarce-lecatism  configuration  determined  in  seep  3  i w  the  stsaa 
configuration  that  was  last  used  in  step  2,  the  iterative  pmesdure 
fcerainates,  and  the  optical  source-1  vocation  configuration  and  den&nd 
allocation  are  approximated  by  the  source-location  configuration  end 
demand  allocation  . orresponding  to  the  current  ^fs  and  B^’g  respectively 
On  the  ether  hand,  if  the  source-location  configuration  has  changed, 
ret  ,ra  to  step  2  with  the  current  values  being  those  just  determined 
In  step  3 

An  extremely  attractive  feature  of  this  iterative  procedure  is  that,  if  the 
objective  function  for  the  general  source-location  problem. 


rd 

>1 


\ij  *kij  +  E£«i 


K 


fki  “ki* 


(46) 


was  calculated  each  time  step  2  was  completed  and  each  tise  step  3  was  completed, 
the  result  would  he  a  non increasing  sequence.  This  characteristic  follows 
immediately  since  in  step  2  the  previous  desand  allocation  is  feasible  whan  the 
new  optimal  demand  allocation  is  determined  and  in  step  3  the  previous  source- 
location  configuration  is  feasible  when  the  new  optimal  configuration  is  determined. 
Thus,  each  time  step  2  is  completed  end  each  tics  step  3  is  completed  an  improved 
feasible  solution  to  the  general  source-location  problem  is  obtained  or  a  feasible 
solution  which  is  at  least  as  good  as  its  predecessor.  This  feature  is  quite 
important  for  vary  lsrjs  problems  in  which  a  oe.  :ch  for  the  optimal  solution  is 
often  economically  impossible. 

it 

Since  there  are  only  (I-fl)  possible  source*»locatioc  cc~*£ jgureci.eas,  the 


iterative  procedure  will  terminate  in  a  finite-  cupbu.  ♦  cf  steps  provided  no 
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cource-location  configuration  is  generated  in  step  3  infinitely  often.  In  step 
2,  the  value  of  (46),  the  objective  function  for  the  general  source- location 
problem,  is  minimised  for  the  current  source-location  configuration.  Thus,  since 
the  sequence  of  values  of  (46)  determined  by  the  iterative  procedure  is  a  non- 
increasing  sequence,  the  only  Hay  in  which  a  source-location  configuration  could 
be  generated  in  step  3  infinitely  often  Is  for  there  to  be  more  than  one  source* 
location  configuration  with  the  same  minimum  value  of  (46)  end  for  looping  or 
cycling  to  occur  among  these  configurations.  However,  such  looping  can  be 
easily  dealt  with.  For  instance,  in  the  examples  to  follow  the  stopping  rule 
in  step  4  was  augmented  sc  that  termination  would  occur  if  the  same  value  of  (46) 
occurred  more  then  four  times ,  Interestingly,  the  iterative  procedure  was  never 
terminated  for  this  reason. 

Since  this  iterative  procedure  yields  only  a  locally  optimal  solution  as 
opposed  to  a  necessarily  globally  optimal  solution,  it  will  generally  be 
worthwhile  to  repeat  the  iterative  procedure  with  different  initial  source- 
locatien  configurations. 

5.  Examples  of  the.  Iterative  Procedure’s  Performance 

In  the  previous  section  several  attractive  characteristics  of  the  iterative 
procedure  were  indicated.  To  further  substantiate  this  attractiveness,  the 
iterative  procedure  was  tried  on  three  sample  problems  involving  real  data.  The 
parameters  anti  problem  characteristics  are  indicated  in  Table  1. 

In  Examples  1  and  2  the  iterative  procedure  was  carried  out  for  each  of  1G0 
initial  source-location  configurations  selected  at  random.  Since  Example  3 
corresponds  to  a  much  larger  sample  problem,  ^netery  considerations  allowed 
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only  33  random  initial  source- location  configurations  to  be  considered. 

the  cost  (46)  associated  with  the  optimal  demand  allocation  for  the  initial 
source- location  configuration  and  the  cost  associated  with  the  optimal  demand 
allocation  for  the  source-lceation  configuration  determined  to  be  approximately 
optimal  by  the  iterative  procedure  were  computed  for  each  of  the  initial  source- 
location  configurations.  As  hoped,  the  costs  associated  with  the  approximate 
solutions  generated  by  the  iterative  procedure  were  much  smaller  than  the  costs 
associated  with  the  optimal  demand  allocations  for  the  random  initial  source- 
location  configurations.  The  magnitude  of  this  improvacnt  is  indicated  in  Table 
2. 

The  effort  involved  in  the  iterative  procedure's  determination  of  an  approximate 
solution  is  reflected  by  the  number  of  source-location  configurations  for  which  the 
corresponding  optimal  demand  allocation  had  to  be  determined.  These  numbers  which 
include  the  initial  source-location  configurations  were  surprlseingly  small. 

As  shown  by  the  summary  of  these  numbers  in  Table  3,  the  iterative  procedure  did 
not  Just  move  gradually  toward  the  near  optimal  solutions  but  rather  leaped  toward 
them. 

.'or  comparative  purposes,  the  sero-one  mixed  integer  linear  progressing 
problem  corresponding  to  the  aesple  problem  in  Example  1  was  also  solved  using  a 
branch-and-bound  algorithm.  This  algorithm  was  initially  progressed  by  Hsstphal 
and  Gately  [26]  who  based  the  algorithm  on  the  branch-and-bound  procedure  described 
by  Davis,  Kendrick,  and  He iceman  [10] .  The  algorithm  was  modified  to  incorporate 
the  improvements  in  the  bounding  procedure  which  were  suggested  by  Tomlin  [23] * 
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The  optimal  solution  has  a  value  of  91  which  is  also  the  yalue  of  the  bsst  approximate 
solution  determined  by  the  iterative  procedure.  The  hranch-and-bound  algorithm 
required  85  minutes  of  IBM  360/65  central  processing  lias  to  determine  the  optical 
solution. 

The  optimal  demand  allocation  for  a  given  source- location  confignration  was 
determined  by  an  out-of-kilter  network-flow  algorithm  described  bT,  Claaen  (27], 

The  total  number  of  IBM  360/65  central  processing  minutes  required  for  ail  of  the 
trials  of  Che  iterative  procedure  in  Example  1,  2,  and  3  was  only  12,  12,  and  20 
respectively. 

6.  Evaluation  of  the  Iterative  Procedure 

The  performance  characteristics  of  the  iterative  procedure  in  the  three 
sample  problems  combined  with  the  attractiveness  of  its  theoretical  sad  practical 
properties  as  described  in  section  4.4  makes  the  iterative  procedure  an  extremely 
attractive  practical  tool  for  determining  near  optimal  solutions  to  general 
source-location  problems. 
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